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ABSTRACT 


Flapping-wing  propulsion  was  studied  experimentally 
through  Laser  Doppler  Velocimetry.  Measurements  were  both 
time-averaged  and  unsteady,  and  were  conducted  on  a  Micro- 
Air  Vehicle  (MAV)  model  developed  at  NPS  by  Professors  Max 
Platzer  and  Kevin  Jones.  The  objective  of  this  work  was  to 
further  understanding  of  the  aerodynamics  of  flapping-wing 
propulsion.  In  specific,  this  study  examined  separation 
control  on  the  leading  fixed  wing  due  to  entrainment  by  the 
trailing  flapping  wings.  Further,  a  study  of  wake  topology 
examined  differences  between  the  optimal  and  off-optimal 
cases.  Experimental  studies  took  place  in  the  NPS  5'  x  5' 
low  speed  wind  tunnel.  The  model  was  supported  on  a  test 
stand  and  LDV  measurements  of  the  flow  field  were  taken. 
Studies  were  made  at  varying  freestream  velocities,  angles 
of  attack,  and  flapping  frequencies.  The  test  stand  was 
instrumented  with  force  balances  to  show  forces  in  both  the 
streamwise  and  vertical  directions. 
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I. 


INTRODUCTION 


A.  OVERVIEW 

The  purpose  of  this  paper  is  to  examine  the  unsteady 
flow  around  a  pair  of  counter-phase  pitching  and  plunging 
airfoils  used  to  power  a  Micro-Air  Vehicle.  At  low  Reynolds 
numbers,  flapping  wings  have  the  potential  to  be  a  more 
efficient  means  of  propulsion  for  small,  lightweight 
aircraft  than  either  conventional  propellers  or  rotary- 
wings.  The  hope  is  that  a  better  understanding  of  the  flow 
field  about  the  flapping  wings  will  allow  for  more 
efficient  design  and  application  of  this  ripening 
technology . 

Experimental  observations  were  made  in  the  NPS  low- 
speed  in-draft  wind  tunnel  with  a  5'  x  5'  test  section,  and 
a  TSI  two-channel  Laser  Doppler  Anemometer. 

B .  BACKGROUND 

The  obviousness  of  effective  propulsion  via  flapping 
wings  was  regarded  as  early  as  Icarus  and  Daedalus,  and 
perhaps  earlier,  though  with  sometime  disastrous  if 
fictional  results.  With  regard  to  a  more  factual  approach, 
any  child  can  recognize  that  a  variety  of  birds  and  insects 
propel  themselves  through  the  air  via  flapping  wings. 
Indeed,  in  the  early  days  of  human  flight,  a  number  of 
attempts  were  made  to  emulate  the  flapping  of  the  birds  as 
a  means  of  flight;  most  meeting  with  singularly  Icaran 
results . 
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Current  wisdom  would  suggest  that  these  early  attempts 
failed  for  any  number  of  reasons  but  especially  the 
structural  weight  of  large  man-carrying  machines  that  must 
also  flap  their  primary  lifting  surfaces.  The  flapping  wing 
MAV  designed  by  Professors  Platzer  and  Jones  and  examined 
herein  addresses  this  problem  in  a  unique  fashion.  First, 
being  a  micro-air  vehicle,  the  machine  is  penalized  neither 
by  the  weight  of  a  pilot  nor  by  a  very  large  structural 
weight.  Even  so,  it  is  not  the  primary  lifting  surface  that 
is  pitched  and  plunged  to  achieve  the  flapping  motion,  but 
rather  two  smaller  trailing  edge  counter-phase  flappers. 
The  purpose  of  this  experiment  was  to  examine 
experimentally  the  flow  field  about  those  two  flapping 
wings  and  the  effect  of  their  flapping  on  the  flow  over  the 
main  wing. 


C .  OBJECTIVES 

The  ostensible  goal  of  this  work  was  to  contribute  to 
the  body  of  knowledge  being  developed  at  NPS  and  elsewhere 
with  regard  to  flapping-wing  propulsion.  In  specific,  the 
primary  goal  was  to  develop  a  time-unsteady  model  of  the 
flow  field  about  the  trailing  edge  pitching  and  plunging 
device  as  developed  by  Professors  Platzer  and  Jones. 
Special  attention  was  paid  to  flow  entrainment  about  the 
leading  fixed  wing  with  respect  to  stall  delay.  Also  of 
interest  were  the  variations  in  wake  patterns  among  off- 
and  on-optimal  cases.  To  this  end,  a  procedure  was 
developed  using  frequency-resolved  LDV  data  to  examine  the 
flow  around  the  full-scale  device  in  a  wind  tunnel 
environment . 
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Ill .  EXPERIMENTAL  APPROACH 


A.  MODELS 

All  models  studied  in  this  report  are  the  product  of 
years  of  analysis  and  experimentation  by  Professor  Kevin 
Jones,  who  is  responsible  for  their  construction  as  well  as 
that  of  the  hardware  used  to  mount  and  instrument  the 
models . 

1 .  Main  Wing  Section 

Though  absent  for  much  of  the  experimental  work,  the 
main  wing  is  crucial  in  generating  the  lift  required  to 
hold  the  operational  MAV  aloft.  Consequently,  the  fixed 
wing  is  added  to  the  wind-tunnel  model  in  order  to  examine 
interaction  between  itself  and  the  trailing  edge  flapping 
wings . 

The  main  wing' s  frame  is  composed  entirely  of  hand- 
carved  balsa  and  carbon  fiber  assembled  with  glue.  The  skin 
is  of  Japanese  tissue  paper.  The  fuselage  houses  the  motor, 
transmission,  battery  and  radio  gear  in  the  operational 
model,  but  only  the  motor  and  transmission  in  the  wind 
tunnel  model . 

The  span  of  the  main  wing  is  slightly  wider  than  that 
of  the  flapping  wings  at  10.5"  with  a  mean  geometric  chord 
of  5.5"  and  absent  both  taper  ratio  and  sweep.  The  main 
wing  also  has  8.5°  of  dihedral. 

The  main  wing' s  airfoil  section  might  ideally  be 
infinitely  thin;  however,  structural  considerations  deny 
this  option  [cf .  Ref  2] . 


3 


2 .  Flapping  Wings 

The  flapping  wings  are  powered  by  a  small  electric 
motor  run  on  an  external  power  supply  for  the  purpose  of 
this  experiment.  The  plunging  motion  of  the  flapping  wings 
is  enforced  by  the  motor's  rotation,  but  the  pitching 
motion  is  aeroelastic  in  nature  and  is  limited  by  the 
stiffness  in  the  flexible  wing  mounts.  The  wings  move  in 
counterphase  to  cancel  out  the  inertial  and  aerodynamic 
effects  of  a  single  wing  moving  up  and  down. 

The  flapping  wings  are  identical  in  construction,  each 
having  a  chord  of  1.5"  and  a  span  of  9.75".  The  leading 
edge  is  a  balsa  dowel.  The  wing's  surface  is  of  microfilm 
interspersed  at  intervals  with  chordwise  battens.  The 
plunging  amplitude,  from  the  centerline  to  top-dead-center, 
is  0.6"  or  0.4  chordlengths ,  with  a  mean  separation  of 
3.2". 


Figure  1 .  Side  View  of  Flapping  Wing  Model 
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3 .  Test  Mount 


The  stand  upon  which  the  model  is  placed  is  isolated 


from 

wind 

tunnel 

vibration 

by  being 

attached  not 

to 

the 

floor 

of 

the  tunnel,  but 

through  a 

covered  hole 

in 

the 

floor 

to 

the  more 

stationary  structure  beneath.  This 

test 

stand 

is 

faired 

into  the 

wind  by 

being  shaped 

like  a 

symmetric  airfoil,  and  is  outfitted  with  sensors  to  resolve 
lift  and  drag  in  the  stability  axes,  i.e.,  in  the 
streamwise  and  vertical  directions.  The  mount  also  allows 
the  model  to  be  rotated  about  the  spanwise  or  "pitch"  axis, 
enabling  the  model  to  be  studied  at  varying  angles  of 
attack . 


Oscilloscope  Demonstrating  Square-Wave 
Signal  with  Lift  and  Drag  Gauges  Beneath 
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Figure  2 . 


B .  TEST  EQUIPMENT 

1 .  Wind  Tunnel 

The  experiments  referred  to  herein  were  each  conducted 
in  the  NPS  low  speed  in-draft  wind  tunnel.  This  tunnel 
ingests  air  from  within  the  building  through  a  15'  square 
intake.  Following  a  9:1  contraction,  the  flow  enters  a  5' 
square  test  section.  Flow  speed  is  modulated  by  a  constant- 
speed,  variable-pitch  electric  fan.  The  fan  is  suspended  on 
rubberized  mounts  to  minimize  vibration  in  the  test 
section.  Further,  to  reduce  swirl  and  turbulence  in  the 
tunnel,  honeycomb  screens  have  been  placed  both  before  and 
after  the  test  section,  though  the  downstream  screens  do 
more  to  protect  the  fans  from  runaway  models  than  to 
straighten  the  flow. 

2 .  Laser  Doppler  Velocimeter 

The  primary  technique  used  to  observe  the  velocities 
in  the  flow  field  about  the  MAV  was  laser  Doppler 
velocimetry.  Specifically,  a  TSI  two-channel  LDV  system 
with  a  single  probe  was  used.  The  streamwise  component  of 
the  flow  was  observed  via  the  blue  channel  at  1.92xl0-5  in 
(488  nm)  wavelength,  and  the  vertical  component  via  the 
green  channel  at  a  wavelength  of  2.03xl0"5  in  (514.5  nm)  . 
Both  colored  beams  were  phase  shifted  by  the  Bragg  cell  at 
a  frequency  of  1  MHz.  For  more  detail  on  the  mechanics  of 
how  an  LDV  systems  measures  flow  velocity,  please  see  Lund, 
Appendix  C . 
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Figure  3.  TSI  Laser 

The  probe  was  fixed  to  a  VP-9000  three-dimensional 
traverse,  which  was  controlled  remotely  via  a  Windows  PC 
running  TSI's  Flow  Information  Display  (FIND)  or  PACE 
software . 
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The  FIND  software,  as  its  name  implies,  also 
interpreted  and  displayed  velocity  data  from  the  LDV 
system.  As  well,  FIND  allowed  the  user  to  tweak  various 
settings  of  his  equipment  in  real  time.  Parameters  under 
the  PC's  control  included  the  frequency  of  the  phase  shift, 
and  the  frequency  band  of  interest.  PACE  is  a  very  similar 
program,  except  that  it  incorporates  the  use  of  the 
Rotating  Machinery  Resolver,  explicated  below. 


Figure  4.  TSI's  PACE  for  Windows  GUI 

Seeding  was  generated  via  a  theatrical  fog  machine 
using  Rosco  non-toxic  smoke  fluid.  This  method  generated 
particles  of  between  lxlO-5"  and  2.5xl0~3"  diameter  [cf.  Ref 
8]  . 
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One  of  the  black  boxes  of  the  LDV  system  was  a 
Rotating  Machinery  Resolver  (RMR) .  This  machine  acquired  a 
lock  on  the  period  of  revolution  of  the  MAV' s  flapping  via 
a  square  wave.  A  shaft  encoder  was  built  to  generate  the 
wave  each  time  the  wings  returned  to  their  centerline 
position.  Due  to  this  rotary  technology,  velocities 
measured  at  different  sectors  in  the  flapping  cycle  could 
be  measured  and  stored  independently  of  each  other.  In  this 
method,  the  flow  field  about  the  flapping  wings  need  not  be 
measured  on  a  time-averaged  basis.  Each  portion  of  the 
cycle  can  be  studied  separately  and  with  relation  to  the 
others . 


9 


c. 


SETUP 


When  fully  assembled,  the  test  apparatus  functionality 
may  be  represented  by  the  flow  chart  in  the  figure  below. 


Figure  5.  Flow  Chart  of  Experimental 

Instrumentation 


An  explanation  of  the  flow  chart  is  in  order. 
Beginning  in  the  upper  left,  the  Voltage  source,  controlled 
by  the  operator,  output  a  DC  current  used  to  power  the  MAV. 
For  the  three  frequencies  used  in  the  experiment,  i.e., 
25.6,  32,  and  38.4Hz,  the  Voltages  needed  by  the  motor  were 


7.9,  9.8,  and  11.8  Volts  respectively.  The  motor  referred 

to  here  is  the  second,  more  powerful  12V  motor  used  for  the 
bulk  of  the  experiment. 

The  rotation  of  the  electric  motor  turned  the 
crankshaft  connected  to  the  flapping  wings.  As  well,  the 
same  shaft  turned  the  optical  encoder,  a  device  built  by 
Prof.  Jones  in  which  a  disk  with  a  single  slot  rotates 
between  an  LED  and  a  photocell.  Each  rotation  of  the  disk 
allowed  a  pulse  of  light  to  pass  through  the  slot  and  to  be 
intercepted  by  the  photocell.  This  cell  then  emitted  a  weak 
electrical  current  that  was  transmitted  to  the  signal 
conditioner . 

The  signal  conditioner,  upon  receiving  the  weak 
electric  pulse  of  the  photocell,  emitted  the  leading  edge 
of  a  short,  5V  square  wave.  This  step  was  necessary  because 
the  RMR,  described  above,  had  very  stringent  requirements 
about  the  signal  it  could  use  to  acquire  a  phase  lock.  It 
needed  a  50Q  TTL  signal. 

The  square  wave  was  also  sent  to  an  oscilloscope,  from 
which  the  operator  read  the  frequency  of  oscillation  of  the 
flapping  wings  and  adjusted  the  voltage  input  accordingly. 
Given  a  constant  voltage,  the  MAV  model  was  quite  good  at 
holding  a  frequency,  to  within  0.25Hz  or  about  0.8%. 

The  newly  conditioned  signal  was  then  sent  to  the  RMR. 

The  RMR  used  the  square  wave  to  acquire  a  phase  locked 

loop.  By  this  method,  the  RMR  can  know  what  part  of  the 

cycle  it  is  in,  even  though  it  only  receives  one  rising- 

edge  input  per  cycle.  That  is  to  say  that  the  machine 

determined  the  time  that  the  MAV  exhausted  in  one  cycle  and 

divided  it  into  360  parts--each  representing  one  degree  of 
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rotation.  When  it  received  the  signal  from  the  Colorburst 
that  a  particle's  velocity  had  been  measured  by  the  LDV, 
the  RMR  made  a  note  of  the  position  where  the  MAV  was  in 
its  cycle  when  the  velocity  was  measured.  The  RMR  then  sent 
a  signal  to  the  IFA  to  flag  that  velocity  with  a  number 
indicating  the  number  of  degrees  elapsed  between  the  rising 
edge  and  the  particle's  passing  through  the  control  volume. 

The  Intelligent  Flow  Analyzer,  model  755,  turned  the 
signal  from  the  Colorburst  into  a  form  usable  by  the  PC  and 
transmitted  the  same  to  the  PC  where  it  was  recorded  by  the 
PACE  software. 

The  Colorburst  was  yet  another  black  box  in  the  LDV 
system.  Colorburst  incorporated  the  phase  shifter  that 
increased  the  velocity  of  particles  in  the  flow  relative  to 
the  lines  of  interference  in  the  probe  volume,  allowing 
measurement  of  negative  velocities.  Its  primary  function, 
however,  was  to  receive  and  digitize  data  received  from  the 
probe  in  the  form  of  f iber-optically  transmitted  pulses  of 
reflected  light. 

The  Argon-Ion  Laser,  of  course,  generated  the  laser 
beams  used  to  create  the  probe  volume  in  the  flow.  Although 
it  is  capable  of  using  much  more  power,  it  was  found  that 
the  best  data  rate  was  achieved  by  using  3.4  BTU/hr  (.99 
Watt)  . 

The  PC  was  very  busy.  Via  the  Windows-based  PACE 
software,  it  was  able  to  control  every  aspect  of  the  LDV 
system.  Using  inputs  from  the  operator,  the  PC  controlled 
the  settings  of  the  Colorburst  and  received  velocity  and 
frequency  data  from  the  IFA.  Following  a  matrix  created  by 
the  operator,  the  PC  was  able  to  determine  when  enough  data 
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had  been  collected  at  a  given  point  in  the  flow  and  to  move 
the  traverse  to  the  next  indicated  point.  The  power  held  by 
the  computer  over  the  LDV  system  enabled  the  experiment  to 
be  almost  completely  automated  after  some  initial  setup. 
However,  it  was  found  to  behoove  the  operator  to  stay 
nearby  and  monitor  the  experiment's  progress,  if  for 
nothing  more  than  to  be  certain  that  the  flow  seeding  was 
properly  adjusted. 


D .  PROCEDURE 

In  the  course  of  designing  the  experiment,  certain 
decisions  appeared,  and  forced  themselves  to  be  dealt  with. 
In  the  interest  of  scientific  purity,  the  following 
paragraphs  will  detail  the  reasons  for  most  of  the  more 
important  decisions. 

Though  it  is  not  known  how  three  dimensional  the  wake 
is  about  the  flapping  wings,  that  three-dimensionality  is 
well  beyond  the  scope  of  this  thesis,  not  only  because  the 
LDV  in  the  low-speed  wind  tunnel  has  only  two  components, 
but  also,  as  cursory  inspection  of  the  appendices  will 
indicate,  even  the  narrow  scope  that  was  applied  to  this 
study  has  generated  massive  amounts  of  data  with  which  to 
be  reckoned  and  written  intelligently  about.  Ergo,  the 
spanwise  location  at  which  to  study  the  flow  field  became  a 
concern.  At  first,  the  centerline  might  seem  a  logical 
place  to  examine  a  wake,  but  the  location  of  the  test  stand 
belied  this  theory.  A  simple  spanwise  survey  of  the  flow 
behind  the  model  generated  the  following  results. 
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Case  1 :  Comparison  of  LDV  Average  to  RMR  Average  at  X=4.5" 


Figure  6.  Spanwise  Survey  of  Turbulence  Intensity 

A  number  of  parameters  were  examined  and  proved 
inconclusive,  including  mean  velocity,  standard  deviation, 
and  seeding  effectiveness.  These  results  tended  to  indicate 
that  the  wake  was  indeed  quite  uniform,  however,  the 
decision  had  to  be  made  and  it  was  finally  based  upon  the 
regularity  of  the  turbulence  intensity  at  1"  to  port  of  the 
centerline . 

That  done,  a  grid  in  the  remaining  two  dimensions  was 
to  be  designed.  It  was  decided  that  the  portions  closer  to 
the  MAV  should  be  spaced  at  k  inch  and  that  the  outer  nodes 
should  be  spaced  at  ^  inch.  Five  columns  were  chosen:  the 
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first  one  chord-length  before  the  flapping  wing' s  leading 
edge,  the  second  immediately  before  the  leading  edge,  the 
third  immediately  aft  of  the  flapping  wing's  trailing  edge, 
the  fourth  one  chord-length  downstream  of  that,  and  yet  a 
fifth  two  more  chord-lengths  downstream.  The  grid  is 
illustrated  in  figure  7 . 


Figure  7  .  Experimental  Grid 

See  the  appendices  for  vector  plots  of  the  velocity  at 
each  node  and  that  velocity' s  variation  throughout  the 
cycle,  normalized  by  the  freestream  velocity. 

The  next  call  to  be  made  involved  which  cases  to 
examine.  Consultation  with  my  advisor  yielded  the  cases 
described  in  Table  1 . 
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Table  1 .  Cases  under  Examination 


Test  Case 

AOA  (deg) 

Vinf  (fps) 

Frequency(Hz) 

Configuration 

1 

15 

9.2 

32 

Wing  attached 

2 

15 

9.2 

0 

Wing  attached 

3 

15 

9.2 

32 

No  fixed  wing 

4 

0 

9.2 

32 

No  fixed  wing 

5 

0 

9.2 

38.4 

No  fixed  wing 

6 

0 

9.2 

25.6 

No  fixed  wing 

One  final,  major  question  remained,  and  this  was  to 
determine  how  much  time  was  spent  and/or  wasted  gathering 
necessary  or  extraneous  data  points.  That  is  to  say, 
exactly  how  many  velocities  need  to  be  measured  by  the  LDV 
before  one  can  say  one  has  indeed  measured  the  velocity.  In 

the  static  case,  this  was  quite  simple  since  1000  points 

can  be  gathered  quite  quickly,  but  for  the  RMR,  each  case 
is  divided  into  72  sections,  and  to  wait  for  72,000  points 
at  each  of  78  nodes  in  5  cases  could  take  an  eternity. 

Consequently,  following  consultation  with  Professors 
Chandrasekhar  and  Hobson,  a  study  was  devised  whereby  the 
relative  merits  of  differing  numbers  of  data  points  could 
be  compared.  The  table  you  see  below  was  the  result,  and  I 
am  thrilled  to  report  that  the  result  was  that  only  3600 
points  (or  50  points  per  5°  of  arc)  were  found  to  be 

necessary . 

Table  2.  Necessary  Number  of  Data  Points 


Avg  pts/5deg 

Min  Points/5deg 

Tl  (%) 

StDev 

50 

25 

3.12 

0.19 

100 

50 

2.93 

0.37 

200 

100 

3.16 

0.28 

400 

200 

3.17 

0.31 
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IV.  RESULTS 


A.  FLOW  FORWARD  OF  THE  FLAPPING  WINGS 

As  you  will  no  doubt  notice  upon  inspection  of  the 
plots,  examination  of  the  flow  ahead  of  the  flapping  wings 
produced  some  rather  fascinating  results.  Witness  Figure  8. 
Here,  the  'x' s  represent  the  air's  horizontal  velocity 
component  immediately  upstream  of  the  flapping  wing' s 
leading  edge,  while  the  'o's  are  one  chord  length  upstream 
of  that.  In  less  than  an  inch  and  a  half,  the  flow  has 
accelerated  from  slightly  above  freestream  speed  to  almost 
half  again  the  same  speed.  The  lack  of  definition  in  the 
flow  over  the  upper  surface  is  due,  as  has  no  doubt  been 
noticed,  to  the  absence  of  points  within  nearly  two  inches 
of  the  model's  centerline.  The  dihedral  of  the  fixed  wing 
prevented  the  LDV' s  probe  volume  from  descending  closer  to 
the  wing's  surface.  This  was  a  problem  for  both  Cases  1  and 
2,  though  for  the  final  four  cases,  the  wing  was  removed 
and  the  probe  allowed  to  descend  significantly  farther. 
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Case  1 :  Flow  over  the  Fixed  Wing 


Figure  8.  Flow  over  the  Fixed  Wing,  Case  1 

By  way  of  contrast.  Figure  9  illustrates  the  same 
model  in  the  identical  flight  condition,  except  that  the 
flapping  wings  are  still  and  resting  in  their  mean 
position.  Note  the  sharp  deceleration  of  the  flow  near  the 
stagnation  point  of  the  flapping  wing's  (now  stationary) 
leading  edge. 
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Case  2:  Flow  over  the  Fixed  Wing 


Figure  9.  Flow  over  the  Fixed  Wing,  Case  2 

Figure  10  is  an  illustration  of  Case  3,  in  which  the 
conditions  are  identical  to  Case  1,  save  for  the  absence  of 
the  fixed  wing.  The  similarity  of  the  plots  indicates  that 
the  presence  or  absence  of  a  fixed  wing  is  no  prerequisite 
for  the  entrainment  of  flow  by  the  trailing  flapping  wings. 
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Case  3:  Flow  forward  of  the  Flapping  Wings 


Figure  10.  Flow  forward  of  the  Flapping  Wings, 

Case  3 

Figure  11  illustrates  how  much  closer  to  the  model's 
centerline  the  probe  was  able  to  approach  in  the  absence  of 
the  fixed  wing.  With  regard  to  the  obvious  asymmetry  above 
and  below,  it  may  be  speculated  that  this,  too,  is  a 
reflection  of  the  mean  position  of  the  flapping  wing  that 
is  simply  missed  by  the  resolution  on  the  upper  surface  and 
was  only  observed  here  by  chance.  It  is  the  author's 
opinion  that  this  merits  further  study. 
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Case  4:  Flow  forward  of  the  Flapping  Wings 


Figure  11.  Flow  forward  of  the  Flapping  Wings, 

Case  4 

Case  5,  in  Figure  12  seems  better  able  to  illustrate 
the  supposed  symmetry  of  the  flow  above  and  below  the 
wings,  while  Case  6,  the  with  the  fastest  flapping  rate  of 
38.4  Hz,  seems  again  to  have  resolved  a  similar  phenomenon, 
perhaps  again  by  chance,  and  definitely  again  on  the  lower 
surface . 
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Case  5:  Flow  forward  of  the  Flapping  Wings 


Figure  12.  Flow  forward  of  the  Flapping  Wings, 

Case  5 
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Case  6:  Flow  forward  of  the  Flapping  Wings 


Figure  13.  Flow  forward  of  the  Flapping  Wings, 

Case  6 


B.  LDV  VERSUS  RMR  MEAN  VELOCITIES 


When  the  LDV  software,  that  is  to  say  FIND®  or  PACE® 


solves 

for 

a  mean  velocity  at 

a  given  locus 

in  the 

flow 

field. 

it 

does  so  by  means 

of  averaging 

all  of 

the 

velocities  it  measures  there.  This  means  that,  if  the  flow 


is  cyclical  and  if  that  cyclical  flow  allows  more  seed 
particles  to  pass  through  the  probe  volume  at  a  certain 
part  of  the  cycle  than  at  other  parts,  the  velocity  will  be 
weighted  toward  that  of  the  heavier  seeded  part  of  the 
cycle . 

Since  the  flow  in  the  MAV' s  wake  is  heavily  cyclical, 

it  provides  a  distinct  opportunity  to  study  this  condition. 
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In  order  to  gather  a  non-weighted  average,  RMR  data 
was  taken.  Now,  the  RMR  will  indeed  collect  more  points 
during  certain  parts  of  the  cycle  than  others,  meaning  that 
certain  folders,  each  representing  a  five  degree  arc  of  the 
cycle,  will  be  fuller  than  others.  To  counter  this  effect, 
each  folder  is  quite  simply  weighted  equally  in  what  is 
referred  to  as  the  RMR  mean  velocity.  This  name  is  meant  to 
distinguish  it  from  the  conventional  LDV  mean  velocity, 
although,  of  course,  both  methods  used  the  same  LDV,  and 
indeed,  the  same  data;  the  only  difference  is  in  the 
weighting . 


Case  1 :  Comparison  of  LDV  Average  to  RMR  Average  at  X=0 


Figure  14.  Comparison  of  LDV  and  RMR  Average 

Velocities,  Case  1  at  X=0" 
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Figure  14  is  fairly  characteristic  of  the  difference 
between  LDV  and  RMR  velocities.  The  two  averages  are 
generally  similar  within  a  few  percent,  but  whether  the  LDV 
average  will  exaggerate  the  high  speeds  or  low  speeds  or 
even  tend  to  wash  them  out  for  the  middle  ground  is  rather 
unpredictable,  since,  as  stated  previously,  it  is  heavily  a 
function  of  seeding. 

C.  WAKE  VARIATION  OVER  THE  CYCLE 

As  one  might  expect,  the  velocity  at  a  given  point 
behind  the  flapping  wing  can  vary  considerably  over  the 
course  of  a  360°  cycle.  The  following  plots  illustrate  some 
of  the  variations  in  both  the  horizontal  direction,  where 
the  mean  is  naturally  positive,  and  the  vertical  direction 
where,  in  most  cases,  it  averages  to  nearly  zero. 
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Case  1 :  Cyclical  Variation  in  Velocity  at  X=  0  (in)and  Z=-1 ,375000e+000. 


Figure  15.  Cyclical  Velocity  Profile  at  X=0  & 

Z=-l . 375 

The  example  chosen  here  for  discussion  and  shown  in 
Figure  15  is  particularly  fascinating  because  it  is  at  a 
point  in  space  immediately  aft  of  the  lower  flapping  wing's 
trailing  edge:  a  point  that  the  trailing  edge  passes  twice 
on  each  cycle,  and  a  particularly  well  defined  one  at  that. 
One  can  clearly  see  that  the  down-stroke  creates  a  peak  in 
the  horizontal  velocity  as  a  jet  of  air  is  forced 
downstream.  As  well,  the  same  passage,  at  about  100°  from 
bottom-dead-center,  generates  two  peaks  in  the  vertical 
velocity  component  as  the  sign  of  the  "V"  velocity  suddenly 
changes.  Bottom-dead-center  is  perhaps  more  appropriate  in 
reference  to  the  upper  flapping  wing,  but  in  any  case  it 
refers  to  the  position  in  which  the  wings  are  as  close 
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together  as  they  can  be,  and  BDC  is  used  here  as  the 
reference  position  of  0°.  It  is  especially  interesting  that 
these  extrema  appear  so  clearly  even  though  the  flapping 
wing  is  not  pointed  directly  downstream.  Indeed,  the  wing 
is  inclined  at  15°  geometric  angle  of  attack. 

D .  VELOCITY  EXTREMES 

In  the  interest  of  developing  not  only  knowledge  about 
the  aerodynamics  of  flapping  wing  propulsion,  but  also 
developing  a  workable  Micro  Air  Vehicle,  one  of  the 
purposes  of  this  study  is  to  examine  wind  velocities  near 
the  model  with  an  eye  to  finding  a  suitable  location  for 
control  surfaces  on  the  flying  model.  The  best  location 
will  have  relatively  consistent  flow.  If  there  is  too 
little  airspeed  the  control  surface  will  be  ineffective, 
too  much  and  the  structural  weight  necessitated  will  ground 
the  lightly  built  MAV.  The  plots  shown  below  represent  the 
minimum,  mean,  and  maximum  velocities  measured  at  each 
location.  For  purposes  of  applicability,  only  the  plots 
from  Case  1  are  shown,  since  they  are  most  similar  to 
actual  flight  conditions.  The  remainder  may  be  viewed  in 
the  appendices. 
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Case  1:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  16.  Minimum,  Mean,  and  Maximum  U  Velocities 

at  X=0 ,  Case  1 
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Case  1 :  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  17.  Minimum,  Mean,  and  Maximum  V  Velocities 

at  X=0 ,  Case  1 

E .  WAKE  DISSIPATION 

Because  a  drag  wake  has  a  lower  kinetic  energy  than 
the  surrounding  flow,  it  can  be  expected  to  be  relatively 
more  stable  than  a  thrusting  wake,  which  has  a  higher 
energy  than  the  surrounding  flow.  This  fact  is  demonstrated 
starkly  in  the  figure  below  where  both  the  thrusting  and 
dragging  wakes  are  extremely  distinct  at  the  X=0  position 
corresponding  to  the  region  immediately  aft  of  the  flapping 
wings'  trailing  edges.  However,  at  the  most  downstream 
points,  only  4.5"  later,  while  the  dragging  wake  a  remains 
strong,  the  thrust  wake  is  already  dissipating  into  the 
freestream.  In  the  current  configuration,  each  flapping 
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wing  generates  its  own  thrusting  wake.  As  these  two  wakes 
move  downstream,  they  spread  apart  and  become  less  distinct 
in  the  process.  In  the  plots  below,  it  will  appear  that  the 
lower  wake  from  Case  1  is  descending  while  the  upper  one 
remains  at  a  fairly  constant  height,  and  even  to  follow  its 
brother  slightly.  The  illusion  is  due  to  the  sharp  angle  of 
attack  at  which  the  model  is  set.  The  wakes  only  move  apart 
relative  to  the  axis  of  the  jet  produced  by  the  flapping. 
In  this  case,  the  jet  is  inclined  at  15°  relative  to  the 
axes  of  the  LDV  system. 


Comparison  of  LDV  Data  between  Flapping  and  Static  Cases  (Cases  1  &2)  at  X=0" 


Figure  18.  Streamwise  Wake  Profile  Comparison  at 

X=0 ,  Cases  1  &  2 
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Comparison  ofLDVData  between  Flapping  and  Static  Cases  (Cases  1  &  2)  atX=1.5" 


Figure  19.  Streamwise  Wake  Profile  Comparison  at 

X=  1.5",  Cases  1  &.  2 
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Comparison  ofLDVData  between  Flapping  and  Static  Cases  (Cases  1  &2)  at  X=4.5" 


Figure  20.  Streamwise  Wake  Profile  Comparison  at 

X=4 . 5" ,  Cases  1  &  2 
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V. 


CONCLUSIONS 


With  respect  to  each  of  the  sub-headings  in  the 
results,  a  number  of  conclusions  may  be  drawn. 

Significant  flow  entrainment  is  apparent  when  the 
flapping  wings  are  active.  Indeed,  examination  of  the  wake 
profile  of  the  static  model  at  15°  angle  of  attack  will 
reveal  that  the  upper  flapping  wing  is  completely  masked  by 
the  wake  region  of  the  fixed  wing.  Ergo,  were  it  not  for 
the  entraining  effects  of  the  flapping  wings,  not  only 
would  the  fixed  wing  be  stalled,  but  also  the  upper 
flapping  wing  would  be  in  such  a  separated  flow  condition 
that  it  would  hardly  be  effective  at  all. 

While  conventional  LDV  velocity  averages  may  be 
excellent  for  static  cases,  it  pales  in  comparison  to  the 
RMR  velocity  average  for  a  cyclical  flow  and  should 
probably  not  be  relied  upon  for  more  than  a  rough  estimate 
in  such  a  case . 

Studies  of  wake  variation  over  the  cycle  reveal  that 
the  wake  velocity  varies  greatly  across  the  regime.  A  more 
focused,  more  highly  resolved  analysis  might  reveal  the 
answers  to  some  of  the  new  questions  posed  by  the  study 
here,  especially  with  respect  to  the  question  about  the 
leading  edge  stagnation  point. 

The  plots  of  velocity  extremes  may  be  used  to  identify 
a  potentially  advantageous  location  for  a  control  surface 
placement  on  the  flying  model.  That  decision  will  be  left 
for  the  master  MAV  designer  and  builder.  Prof.  Kevin  Jones, 
for  consideration. 
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The  thrust  wake's  spreading,  as  predicted  by  potential 
flow  theory,  is  apparent  in  the  wake  profile  plots; 
however,  over  the  limited  distance  considered,  it  is 
difficult  if  not  impossible  to  draw  conclusions  about  the 
thrust  wake's  dissipation  compared  to  that  of  the  drag 
wake.  In  this  respect,  the  physical  dimensions  of  the 
traverse  mechanism  upon  which  the  LDV  probe  rests  impose  a 
limitation . 
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APPENDIX  A.  MATLAB  CODES 


A.  DATA  FORMATTING 

This  code  was  used  to  get  the  .  rmr  files  into  a  form 
usable  by  MATLAB.  Each  file  was  first  imported  as  a  matrix 
called  "data,"  and  then  the  code  was  run. 

%load  data  into  matrix,  omitting  empty  W  columns 
combo=zeros (360, 9) ; 

combo ( : , 1 ) =data ( : , 2 ) ;  %encoder  position 

combo ( : , 2 ) =data ( : , 3 ) ;  %U  pts 

combo ( : , 3 ) =data ( : , 4 ) ;  %V  pts 

combo ( : , 4) =3 . 28084*data ( : , 6) ;  %  U  mean  (fps) 

combo ( : , 5) =3 . 28084*data ( : , 7) ;  %  V  mean  (fps) 

combo ( : , 6 ) =data ( : , 9 ) ;  %  U  std 

combo ( : , 7 ) =data ( : , 10 ) ;  %  V  std 

combo ( : , 8) =data ( : , 12) ;  %  U  TI% 

combo ( : , 9) =data ( : , 13) ;  %  V  TI  % 

%Compress  360  rows  into  72  rows 
mat=zeros (72, 9) ; 
for  1=2:12; 

mat (i, : ) =sum (combo ( [ (5*i-6)  :  (5*i-2)],:)); 

end 

mat (1, : ) =sum (combo ([359,360,1,2,3],:)); 
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mat (:,!)  = [0:5:355] 


save  cl_.xls  mat  -ascii  -tabs 
clear 

B .  PLOT  GENERATION 

This  code  was  then  run  to  generate  the  desired 
plots . 

clear  all 

%initialize  memory 
a=zeros (5616, 9) ; 

%define  indices 
s=zeros (72,1) ; 
e=s; 

for  i=l : 7  8 

s (i) =72* (i-1) +1; 
e (i) =72*i; 

end 

%load  data  files  into  the  matrix 
i=l; 

a  (s  (i)  :e(i) , : ) =load  (' cl_l.xls') ; 
i=i+l ; 
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a (s (i) :e(i) , : ) =load  ('cl_2.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_3.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  (' cl_4.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_5.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_6.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_7.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_8.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_9.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_10.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_ll.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_12.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_13.xls') ; 
i=i  +  l ; 


37 


a (s (i) :e(i) , : ) =load  ('cl_14.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_15.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_16.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_17.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_18.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_19.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_20.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_21.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_22.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_23.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_24.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_25.xls') ; 
i=i  +  l ; 
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a (s (i) :e(i) , : ) =load  ('cl_26.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_27.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_28.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_29.xls') ; 
i=i  +  l ; 

a (s (i) :e(i) , : ) =load  ('cl_30.xls') ; 
i=i  +  l  ; 

a (s (i) :e (i) , : ) =load  ('cl_31.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_32.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_33.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_34.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_35.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_36.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_37.xls') ; 
i=i  +  l ; 
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a (s (i) :e(i) , : ) =load  ('cl_38.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_39.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_40.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_41.xls') ; 
i=i  +  l ; 

a (s (i) :e(i) , : ) =load  ('cl_42.xls') ; 
i=i  +  l  ; 

a (s (i) :e (i) , : ) =load  ('cl_43.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_44.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_45.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_46.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_47.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_48.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_49.xls') ; 
i=i  +  l ; 
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a (s (i) :e(i) , : ) =load  ('cl_50.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_51.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_52.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_53.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_54.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_55.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_56.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_57.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_58.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_59.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_60.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_61.xls') ; 
i=i  +  l ; 
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a (s (i) :e(i) , : ) =load  ('cl_62.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_63.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i), : ) =load  ('cl_64.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_65.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_66.xls') ; 
i=i  +  l  ; 

a (s (i) :e(i) , : ) =load  ('cl_67.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_68.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_69.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_70.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_71.xls') ; 
i=i  +  l ; 

a (s (i) :e (i) , : ) =load  ('cl_72.xls') ; 
i=i  +  l ; 

a (s (i) : e (i) , : ) =load  ('cl_73.xls') ; 
i=i  +  l ; 
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a (s (i) :e(i) , : ) =load  ('cl_74.xls') ; 
i=i  +  l ; 

a (s (i) :e(i) , : ) =load  ('cl_75.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_76.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_77.xls') ; 
i=i  +  l ; 

a (s (i) :e(i), : ) =load  ('cl_78.xls') ; 
i=i  +  l ; 

ldv=xlsread ( ' Casel . xls ' )  ; 

%inspect  matrices  to  be  sure  they  are  properly 
%built  by  being  certain  that  columns  are  not  repeated 

%Compare  matrix  to  itself 
for  i=  ( 1 : 78 ) 

if  a  (s  (i)  ,4)==a(s  (i)  ,5) 

disp  ('rebuild  matrix  number') 
disp  (i) 

end 

end 

%Compare  matrices  to  each  other 
for  i= (2 : 78 ) 
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if  a  (s  (i)  ,4)==a(s  (i-1)  ,4) 


disp  ('rebuild  matrix  number',  i) 
elseif  a (s (i) , 5) ==a (s (i-1) , 5) 

disp  ('rebuild  matrix  number',  i) 

end 

end 


%Ensure  minimum  number  of  points  in  each  file 
minupts=min (a (1 : 5616, 2) ) 
minvpts=min (a(l:5616,3)  ) 

%Normalize  data 
Vinf=9.1864  %fps 
a ( : , 4 ) =a ( : , 4 ) /Vinf ;  %  U  mean 
a ( : , 5) =a ( : , 5) /Vinf ;  %  V  mean 


%plots 

zcoords= [3.875  ;3.375;  2.875  ;2.375;  2.125  ; 

; 1.625  ; 1.375;  1.125  ; 0.875  ; 0.625;  0.375  ; 0.125; 

-0.125;  -0.375;  -0.625  ;-0.875  ;-1.125;  -1.375;  - 
; -1.875  ; - 2 . 1 2 5 ;  -2.375;  -2.875  ;-3.375  ; -3.875] ; 


1 . 875 

1 . 625 
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%non-point  weighted  average 
for  i=l:78 

rmr (i,  : ) =mean (a  (s  (i)  :  (e  (i)  )  ,  :  )  )  ; 

Uavg (i, 1) =mean (a (s  (i)  :  (e(i) ) ,4) ) ; 

Vavg (i, 1) =mean (a (s (i) : (e(i) ) ,5) ) ; 

end 

fig=l; 

figure (fig) 

plot (Uavg (1:26,1) , z coords , ' bo- ' ) 
hold  on 

plot (Uavg (27:52,1) , z coords , ' rx- ' ) 
hold  on 

plot (Uavg (53:78,1) , z coords, ' kd- ' ) 
title  ('Case  1:  Wake  Profile,  RMR') 
ylabel (' Vertical  Position  "Z"  (in) ') 
xlabel ( 'Non-dimensional  Velocity' ) 
legend 

legend ( ' x=0 ' , ' x=l . 5 " ' , ' x=4 . 5 " ' ) 


%Comparison  to  point  weighted  average 
f ig=f ig+1 ; 
figure (fig) 


45 


plot(ldv(l:26,5) /Vinf, z coords, 'bo- ' , 


Uavg (1:26) , z coords , ' rx- ' ) 
legend ( ' ldv ' , ' rmr ' ) 

title ('Case  1:  Comparison  of  LDV  Average  to 
Average  at  X=0 ' ) 

ylabel (' Vertical  Position  "Z"  (in) ') 
xlabel ( 'Non-dimensional  Velocity' ) 

f ig=f ig+1 ; 
figure (fig) 

plot (ldv (27:52,5) /Vinf, z coords, 'bo- ' , 

Uavg (27:52) , z coords , ' rx- ' ) 

legend ( ' ldv ' , ' rmr ' ) 

title ('Case  1:  Comparison  of  LDV  Average  to 

Average  at  X=1.5"') 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity' ) 

f ig=f ig+1 ; 
figure (fig) 

plot (ldv (53:78,5) /Vinf, z coords, 'bo- ' , 

Uavg (53:78) , z coords, ' rx- ' ) 

legend ( ' ldv ' , ' rmr ' ) 

title ('Case  1:  Comparison  of  LDV  Average  to 

Average  at  X=4.5"') 


RMR 


RMR 


RMR 


46 


ylabel (' Vertical  Position  "Z"  (in) ') 
xlabel ( 'Non-dimensional  Velocity' ) 


Umin=zeros (78, 1)  ; 

Vmin=zeros (78, 1)  ; 
for  i=l : 7  8 

Umin (i, 1) =min (a (s (i) :e(i) ,4) ) ; 

Vmin (i, 1) =min (a (s (i)  :e(i) ,5) )  ; 

Umax ( i , 1 ) =max ( a ( s ( i )  :e(i)  ,4)  )  ; 

Vmax ( i , 1 ) =max ( a ( s ( i ) :e(i) ,5) ) ; 
end 

f ig=f ig+1 ; 
figure (fig) 

plot (Umin (1:26) , z coords , ' bo- ' ,  Uavg (1:26) , z coords , ' rx- 
' , Umax (1:26) , zcoords , ' kd- ' ) 

legend ('U  min','U  mean',  'U  max',0) 

title ('Case  1:  Comparison  of  Minimum,  Mean,  and 

Maximum  U  Velocities  for  each  Point  at  X=0 ' ) 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity' ) 

f ig=f ig+1 ; 
figure (fig) 
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plot (Umin (27:52) , z coords , ' bo- ' , 


Uavg (27:52) , z coords , ' rx- ' , Umax (27:52) , z coords , ' kd- ' ) 
legend ('U  min','U  mean',  'U  max',0) 

title ('Case  1:  Comparison  of  Minimum,  Mean,  and 

Maximum  U  Velocities  for  each  Point  at  X=1.5"') 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity' ) 

f ig=f ig+1 ; 
figure (fig) 

plot (Umin (53:78) , z coords, 'bo- ' , 

Uavg (53:78) , z coords , ' rx- ' , Umax (53:78) , z coords , ' kd- ' ) 

legend ('U  min','U  mean',  'U  max',0) 

title ('Case  1:  Comparison  of  Minimum,  Mean,  and 

Maximum  U  Velocities  for  each  Point  at  X=4.5"') 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity' ) 


f ig=f ig+1 ; 
figure (fig) 

plot (Vmin (1:26) , z coords , ' bo- ' ,  Vavg (1:26) , z coords , ' rx- 
' , Vmax (1:26) , zcoords , ' kd- ' ) 

legend('V  min','V  mean',  'V  max',0) 

title ('Case  1:  Comparison  of  Minimum,  Mean,  and 

Maximum  V  Velocities  for  each  Point  at  X=0 ' ) 
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ylabel (' Vertical  Position  "Z"  (in) ') 
xlabel ( 'Non-dimensional  Velocity' ) 

f ig=f ig+1 ; 
figure (fig) 

plot (Vmin (27:52) , z coords , ' bo- ' , 

Vavg (27:52) , z coords , ' rx- ' , Vmax (27:52) , z coords , ' kd- ' ) 

legend('V  min','V  mean',  'V  max',0) 

title  ('Case  1:  Comparison  of  Minimum,  Mean, 

Maximum  V  Velocities  for  each  Point  at  X=1.5"') 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity'  ) 

f ig=f ig+1 ; 
figure (fig) 

plot (Vmin (53:78) , z coords, 'bo- ' , 

Vavg (53:78) , z coords, ' rx- ' , Vmax (53:78) , z coords, ' kd- ' ) 

legend('V  min','V  mean',  'V  max',0) 

title ('Case  1:  Comparison  of  Minimum,  Mean, 

Maximum  V  Velocities  for  each  Point  at  X=4.5"') 

ylabel (' Vertical  Position  "Z"  (in) ') 

xlabel ( 'Non-dimensional  Velocity' ) 


zpos= [ zcoords ; zcoords ; zcoords ]  ; 


and 


and 
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for  i=l : 7  8 


%define  x  location 
xpos=l . 5 ; 
if  i<27 

xpos=0 ; 
elseif  i>52 
xpos=4 . 5 ; 

end 

f ig=f ig+1 ; 
figure (fig) ; 

plot (a (s  (i)  :e (i) ,4) , a  (s (i)  :e(i)  ,1)  ,  'bo-')  %U 
hold  on 

plot  (a  (s  (i)  :e(i)  ,5)  ,  a  (s  (i)  :e(i)  ,1)  ,  '  rx-  '  )  %V 
legend ( 'U' , 'V' , 0) 
z=zpos (i, 1) ; 

string=sprint f ( ' Case  1:  Cyclical  Variation  in  Velocity 
at  X=  %d  (in) and  Z=%d.  %d  (in)%d.',  xpos,  z); 
title (string) 

ylabel (' Degrees  from  Bottom  Dead  Center') 

xlabel ( 'Non-dimensional  Velocity' ) 

end 
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%Save  the  plots  to  file 
for  i=l : fig 
if  i>10 

fname=sprintf ( ' Clpt%d . jpg%d. ' , i-10 )  ; 
elseif  i==l 

fname=sprintf ( ' Clwakpro . jpg ' ) ; 
elseif  i==2 

fname=sprintf ( ' CllrO . jpg ' ) ; 
elseif  i==3 

fname=sprintf ( ' Cllrl5 . jpg ' ) ; 
elseif  i==4 

f name=sprint f ( ' Cllr45 . jpg ' ) ; 
elseif  i==5 

fname=sprintf ( ' ClUmmmO ' ) ; 
elseif  i==6 

f name=sprint f ( ' ClUmmml5 . jpg ' ) ; 
elseif  i==7 

f name=sprint f ( ' ClUmmm45 . jpg ' ) ; 
elseif  i==8 

f name=sprint f ( ' CIVmmmO . jpg ' ) ; 
elseif  i==9 

f name=sprint f ( ' ClVmmml5 . jpg ' ) ; 
elseif  i==10 

fname=sprintf ( ' ClVmmm45 . jpg ' ) ; 


51 


end 


print (figure (i) , ' -djpeg90 ' , ' -rlOO ' , fname) 

end 

%Save  data  to  files 
for  i=l:78 

ptmat=a ( s  ( i )  : e  ( i ) ,  :  )  ; 
fname=sprintf ( ' Clpt%d .xls%d. ',i) ; 
save  (fname,  'ptmat -ascii -tabs ' ) 

end 

disp ( ' Done . ' ) 
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APPENDIX  B.  PLOTS  FROM  CASE  1 


Case  1:  Wake  Profile,  RMR 


Figure  21.  Wake  Profiles  for  Case  1 


Case  1 :  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  22.  Minimum,  Mean,  and  Maximum  U  Velocities 

for  X=0,  Case  1 
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Case  1 :  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  23.  Minimum,  Mean,  and  Maximum  V  Velocities 

for  X=0,  Case  1 


Case  1:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=1.5" 


Figure  24  . 


Minimum, 

for 


Mean,  and  Maximum  U  Velocities 
X=1.5",  Case  1 
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Case  1:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  atX=1.5" 


Figure  25.  Minimum,  Mean,  and  Maximum  V  Velocities 

for  X=1.5",  Case  1 


Case  1:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=4.5" 


Figure  26.  Minimum,  Mean,  and  Maximum  U  Velocities 

for  X=4.5",  Case  1 
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Case  1:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=4.5" 


Figure  27.  Minimum,  Mean,  and  Maximum  V  Velocities 

for  X=4.5",  Case  1 
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APPENDIX  C.  DATA  FROM  CASE  2 


Case  2:  Steady  State  Drag  at  X=  0(in) 


Figure  30.  Drag  Wake  Profile  at  X=0 


Case  2:  Steady  State  Drag  at  X=  2.500000e-001(in) 


Figure  31.  Drag  Wake  Profile  at  X=0.25" 
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Case  2:  Steady  State  Drag  at  X=  5.000000e-001  (in) 


Figure  32.  Drag  Wake  Profile  at  X=0.5" 


Case  2:  Steady  State  Drag  at  X=  7.500000e-001(in) 


Figure  33.  Drag  Wake  Profile  at  X=0.75" 
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Case  2:  Steady  State  Drag  at  X=  1  (in) 


Figure  34.  Drag  Wake  Profile  at  X=l" 


Case  2:  Steady  State  Drag  at  X=  1 . 250000 e +000 (in) 


Figure  35. 


Drag  Wake 


Profile  at  X=1.25" 
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Case  2:  Steady  State  Drag  at  X=  1 ,500000e-K)00(in) 


Figure  36.  Drag  Wake  Profile  at  X=1.5" 


Case  2:  Steady  State  Drag  at  X=  2(in) 


Figure  37.  Drag  Wake  Profile  at  X=2" 
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Case  2:  Steady  State  Drag  at  X=  2.500000e-K)00(in) 


Figure  38.  Drag  Wake  Profile  at  X=2.5" 


Case  2:  Steady  State  Drag  at  X=  3(in) 


Figure  39.  Drag  Wake  Profile  at  X=3" 
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Case  2:  Steady  State  Drag  at  X=  4.500000e-K)00(in) 


Figure  4  0 . 


Drag  Wake 


Profile  at  X=4.5" 
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APPENDIX  D.  PLOTS  FROM  CASE  3 


Case  3:  Wake  Profile,  RMR 


Figure  41.  Wake  Profiles  for  Case  3 


Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  42  . 


Minimum,  Mean,  and  Maximum  U  Velocities 
at  X=0 ,  Case  3 
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Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  43. 


Minimum,  Mean, 
at  X=0, 


and  Maximum  V  Velocities 
Case  3 


Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=1.5" 


Figure  44  . 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=1.5",  Case  3 
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Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=1.5" 


Figure  45. 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=1.5",  Case  3 


Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=4.5" 


Figure  46. 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=4 . 5" ,  Case  3 
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Case  3:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=4.5" 


Figure  47  . 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=4.5",  Case  3 
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APPENDIX  E.  PLOTS  FROM  CASE  4 


Case  4:  Wake  Profile,  RMR 


Figure  50.  Wake  Profiles  for  Case  4 


Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  51.  Minimum,  Mean,  and  Maximum  U  Velocities 

at  X=0 ,  Case  4 
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Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  52 . 


Minimum,  Mean, 
at  X=0, 


and  Maximum  V  Velocities 
Case  4 


Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=1.5" 


Figure  53. 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X= 1.5",  Case  4 
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Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=1.5" 


Figure  54  . 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=1.5",  Case  4 


Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=4.5" 


Figure  55. 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=4 .5",  Case  4 
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Case  4:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=4.5" 


Figure  56. 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=4.5",  Case  4 
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3 


4>=  315 

jure  58  . 


4>=  330  <|)=  345 

Case  4  Wake  Profiles,  180-34^ 


APPENDIX  F.  PLOTS  FROM  CASE  5 


Case  5:  Wake  Profile,  RMR 


Figure  59.  Wake  Profiles  for  Case  5 


Case  5:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  60.  Minimum,  Mean,  and  Maximum  U  Velocities 

at  X=0 ,  Case  5 
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Case  5:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  61 . 


Minimum,  Mean, 
at  X=0, 


and  Maximum  V  Velocities 
Case  5 


Case  5:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=1.5" 


Figure  62 . 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=1.5",  Case  5 
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Case  5:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=1.5" 


Figure  63.  Minimum,  Mean,  and  Maximum  V  Velocities 

at  X=1.5",  Case  5 


Case  5:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=4.5" 


Figure  64 . 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=4 . 5" ,  Case  5 
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Figure  65. 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=4 . 5" ,  Case  5 
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(j)=22  5  (>=24  0  (>=255 


<>=27  0  <>=285  <>=30  0 


<>=330  <>=34  5 


Case  5  Wake  Profiles,  180-345 
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APPENDIX  G.  PLOTS  FROM  CASE  6 


Case  6:  Wake  Profile,  RMR 


Figure  68.  Wake  Profiles  for  Case  6 


Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=0 


Figure  69.  Minimum,  Mean,  and  Maximum  U  Velocities 

at  X—  0 ,  Case  6 
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Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=0 


Figure  70.  Minimum,  Mean,  and  Maximum  V  Velocities 

at  X=0 ,  Case  6 


Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=1.5" 


Figure  7 1 . 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=1.5",  Case  6 
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Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=1.5" 


Figure  72  . 


Minimum, 

at 


Mean,  and  Maximum  V  Velocities 
X=1 . 5" ,  Case  6 


Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  U  Velocities  for  each  Point  at  X=4.5" 


Figure  73. 


Minimum, 

at 


Mean,  and  Maximum  U  Velocities 
X=4 . 5" ,  Case  6 
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Case  6:  Comparison  of  Minimum,  Mean,  and  Maximum  V  Velocities  for  each  Point  at  X=4.5" 


Figure  7 4 . 


Minimum, 

at 


Mean,  and  Maximum 
X=4 . 5" ,  Case  6 


Velocities 
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